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ACCOMPLISHMENTS 
Hydrostatic  Foil  Bearings 

Hydrostatic  planar  foil  bearings  have  been  analyzed  using  lubrication 
theory.  The  results  have  been  summarized  in  [1]. 

The  effects  of  fluid  inertia  in  this  type  of  bearing  have  been  analyzed 
and  shown  to  be  of  essential  importance.  The  results  are  given  in  a  technical 
report  included  herewith  in  part  2  of  this  report.  (See  Item  [2]  above). 

Stick  Slip  Vibrations  of  Foils 

A  model  has  been  developed  which  allows  simulation  of  longitudinal 
vibrations  in  foils  under  the  influence  of  Coulomb  friction  with  variable  and 
discontinuous  friction  coefficient.  Self-excited  vibrations  and  stick-slip  phenomena 


have  been  described.  The  results  are  summarized  in  item  [3]  above.  Coupling 
with  transverse  vibrations  has  been  studied  but  not  completed. 

Foil  Response  to  Bumps 

Analysis  has  been  matte  of  the  effect  of  a  rotating  bump  (representing  a 
magnetic  head)  on  a  foil  floated  on  a  cylindrical  drum.  The  effect  of  several 
parameters  have  been  studied  and  in  particular  the  ratio  of  bump  velocity  to  the 
speed  of  propagation  of  transverse  disturbances  in  the  foil.  The  results  are  given 
in  item  [4]  above.  Additional  analyses  and  experiments  including  the  effects  of 
cutouts  have  been  performed  but  have  not  reached  publication. 
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FOIL  BEARING  INVESTIGATION 
Technical  Report:  Anelysis  of  Hydrostatic  Foil  Bearings 
Including  the  Effects  of  Fluid  Inertia 

INTRODUCTION 


I 
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This  report  explains  a  peculiarity  found  in  earlier  analyses  of  externally 
pressurized  foil  bearings.  Ref.  [1]  summarizes  some  of  the  preceding  studies  [2-5]. 

In  the  above  mentioned  paper,  the  problem  of  a  planar  hydrostatic 
foil  bearing  was  analyzed  on  the  basis  of  lubrication  theory.  Among  the  results, 
the  paper  presented  charts  of  foil  bearing  inlet  pressure  vs.  flow  rate.  An  interesting 
feature  of  these  charts  (Figure  5f  is  that  for  the  same  inlet  pressure  a  dual  value 
of  flow  rate  is  predicted.  This  result  is  not  physically  impossible  (since  the  restrictor 
characteristics  would  generally  select  the  proper  unique  physical  solution),  but 
nevertheless,  it  is  unexpected.  The  motivation  for  the  present  work  was  the  hypo¬ 
thesis  that  the  source  of  the  phenomenon  is  fluid  inertia  effects  which  were  not 
included  in  the  analysis.  Indeed,  as  the  results  of  the  present  study  show,  inertia 
effects  alter  drastically  the  behavior  at  higher  flow  rates  and  resolve  the  puzzle 
associated  with  the  earlier  results. 

The  model  under  consideration  is  shown  in  Figure  1.  It  describes  a 
planar  flow  through  a  restrictor  feeding  into  a  cylindrical  region  covered  by  a 
tensioned  foil.  The  pressure  vs  flow  characteristics  of  the  configuration  comprise 
the  objective  of  the  analysis. 

Basic  Equations 

The  starting  point  of  the  analysis  is  the  compressible  Navier-Stokes 
equations  in  polar  coorindates: 


•See  Fig.  2  of  the  present  report,  for  A  *  0. 
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The  continuity  equation  is 


(3) 

apop1)  ft(p*f) 
ftr  +  d* 


The  perfectly  flexible  foil  equation  is 
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The  isothermal  equation  of  state  is 
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The  curvature  in  polar  coordinates  is 


(5) 


(6) 


and  the  clearance 
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The  pressure  boundary  conditions  are  that  the  pressure  on  the  foil 
becomes  negligible  sufficiently  far  from  the  spindle.  No  velocity  boundary  condi¬ 
tions  are  prescribed  as  explained  in  [6]. . 

Simplification 

We  introduce  a  dimensionless  small  parameter  e  to  be  defined  later 
whose  purpose  is  to  indicate  how  open  the  bearing  is  to  flow.  Qualitatively,  it 
is  expected  to  be  directly  related  to  M,  the  mass  flow  rate  and  inversely  related 
to  the  tension.  In  addition  we  define 


T/r, 


compressibility  parameter 


and  the  dimensionless  variables 
*  OZ.  ^ 


it  s 


£L 


"ftt r-) 


<rr=~^~ 
v  -  >ft 


where  k,  m,  n  are  constants  to  be  determined.  »te  icwnn 


(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 


Figure  1  Cylindrical  Hydrostatic  Foil  Bearing 
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Substituting  the  notation  (8)  -  (14)  in  Eqns.  (1)  -  (7)  two  dimensionless 
parameters  arise 


r  _  r 
t  n.-r 


Foil  Bearing  Parameter 


(15) 
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-Zrt 


p  rT 
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Inertia  Parameter 


(16) 


In  accordance  with  conventional  foil  bearing  theory  [7] ,  the  following 
requirements  are  imposed  on  the  equations  for  the  purpose  of  the  determination 
of  k,  m,  n  as  the  limit  -*a  is  pursued. 


1.  The  expansion  of  equation  (6)  has  to  retain  the  capability  of  satis¬ 
fying  the  condition  that  the  tape  becomes  flat  at  infinity. 


2.  The  purpose  of  this  work  is  to  study  cases  in  which  inertia  effects 
are  significant.  It  will  be  required,  therefore,  that  at  least  some  inertia  term  will 
be  of  the  order  of  some  viscous  term. 


3.  Finally,  it  is  required  that  the  lateral  and  longitudinal  terms  in  the 
continuity  equation  will  be  of  the  same  order. 


These  three  requirements  results  in  the  conditions 

;  kxl.  ”7' 


With  these  results,  the  equations  can  be  shown  to  degenerate,  to  the  zeroth  approx¬ 
imation,  into  Prandtl's  boundary-layer  equations. 
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We  assume  that  the  velocity  profile  is  parabolic.  The  profile  which 


satisifes  the  boundary  conditions 

[iso  'Vj-o 

(21) 

’T]=H 

(22) 

(23) 

where  f,  a  measure  of  local  velocity,  is  to  be  determined. 

The  three  unknown  functions  w  will  be  determined  from  the 
momentum  and  continuity  integrals  between  'j^o  and  Vj—ff  and  from  the 
elastic  equation.  ' 


The  momentum  integral  is 

Noting  that  the  integral  of  Eqn.  (3)  with  respect  to  dr  yields  the  mass  flow  rate 
M,  the  continuity  integral  becomes 

[ 'cfrr)£H  _ 

&  ^  «» 

Eliminating  j-  and  7f  from  Eqns.  (20),  (24),  (25)  the  formulation  is 
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together  with  the  boundary  conditions: 


(26) 
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outside  the  lubrication  zone; 
origin  of  ^  is  at  point  of  tengency 


(27) 

(28) 


and 


3 


(atthe  supply  line  boundary;  ^-^) 

<7 


(29) 


(30) 


= T 

The  solution  is  expected  to  depend  on  I,  C,  ^  .  i.e.,  the  problem 

depends  on  four  parameters.  In  the  incompressible  (C-£><x>)  inertialess  (I  *  0)  limit, 
the  formulation  degenerates  to  that  of  [1], 


It  is  useful  to  present  the  results  as  dimensionless  pressure  vs  flow  rate. 
Therefore,  an  alternative  parameter  will  be  defined  which  eliminates  flow  rate 


(31) 


Hence,  the  solution  will  be  of  the  form 


(32) 


(33) 
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RESULTS 


Due  to  extraneous  constraints  the  results  described  here  are  limited  to 
some  particular  cases.  A  desirable  extensive  parametric  investigation  has  not  been 
made.  Nevertheless,  the  main  point  of  the  study,  i.e.,  the  essential  importance  of 
inertia  effects  is  well  demonstrated  in  Figure  2.  It  is  seen  that  for  low  flow  rate 
the  solution  tends  to  be  identical  with  that  of  Ref  [1].  For  higher  flow  rates  the 
behavior  is  completely  different,  showing  that  increased  groove  pressure  is  required, 
as  one  may  expect,  in  order  to  overcome  additional  inertia  effect. 
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APPENDIX 


The  solution  procedure  is  described  below.  Let  equations  (26M28)  be 

A  a 

reformulated  in  terms  of  alternative  variables  H(t),  in  order  to  permit  solution 
as  an  initial  value  problem:  (Later  the  solution  for  H  will  be  transformed  back  into 
H  ).  After  studying  the  procedure  the  reader  may  verify  that  the  following 
formulation  may  indeed  be  recast  into  the  original  coordinates  using  transformation 
formula  given  below.) 


m 

H  - 


(  -  otA« 


-a 


/>  If 


(33) 


Subject  to  the  initial  conditions,  that  at 
A  * 

A  /  Gt»  / 

H  =  rtf 

A" 

*  -  J 
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A  AU 


(34) 

(35) 

(36) 


where  the  parameters  o/*?  ol  are  inputs.  Integrating  numer- 

ically  from  ,  then  depending  on  the  selection  of  inputs,  the  initially  positve 

y  A 

H  may  become  positively  unbounded  for  large  ^  or  it  may  decrease  to  zero  causing 

the  denominator  of  (33)  to  vanish.  Only  the  first  event  is  of  physical  interest.  In 
A  A 

this  case,  H"  tends  to  reach  an  asymptotic  value  denoted  by  H"  .  Once  the 
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numerical  value  of  this  quantity  has  reached  sufficiently  close  to  its  asymptotic 

value,  the  solution  may  be  terminated  (at  a  somewhat  arbitrary  point)  with  the 
A  \  r\ 

known  values  of  H  ,  H  W  retained.  Now  the  physical  parameters  of  the  problem 


h=f<  £ 


Vi  VjT 


(37) 


1=  (*'l 


Vr  ,x 
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(38) 


where 

*  * 
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While  the  formulation  (26)  -  (30)  depends  on  four  parameters  shown  in  Eqn.  (31), 
the  equivalent  formulation  (33)  -  (36)  depends  on  six  inputs.  Thus,  two  of  the 
six  parameters  are  expected  to  be  arbitrary,  and  to  lead  to  redundant  solutions. 


The  iterative  approach  selected  is  based  on  consideration  of  the  following 
limiting  cases.  By  comparing  Eqns.  (33)  and  (26)  one  may  see  that  with  ■  1, 
the  limit  i/j’O  leads  to  the  incompressible  case,  and  the  limit  •  0  leads  to  the 
inertia-less  case.  Finally,  the  case  of  a  groove  located  at  the  center  of  the  lubrication 
zone  ( ■  ^  )  requires  (Eqn.  29)  H0  *  0.  Thus,  one  may  expect  that  to  a  first 
approximation  He  \  ^  control  £1  /££,.  ^ 

respectively.  To  determine  the  fourth  parameter  e.g.,  £  (or  flow  rate)  we  choose 

A  C.  A  ti 

Hg".  Thus  we  arbitrarily  fix  <xR*l  and  Hg  *  1,  then  for  each  selected  value  of  Hg, 

we  iterate  on  H  ',  o(  o(  to  generate  a  solution  with  desired  values  of  J-  A 
9  *7  *>  0  /  C-  > 

and  some  point  on  the  ^  -  ^  curve. 
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